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Reliability  oaoers,  technical  reports,  and  textbooks  almost  univer¬ 
sally  document  the  derivation  of  availability  equations  for  two  unit  full- 
on  parallel  redundant  systems  with  the  assumption  of  "multiple  ^epai^" 
crews  or  facilities.  Multiple  repair  is  the  case  where  there  are  as  many 
repair  crews  available  as  there  are  possible  units  to  fail.  Some  reli¬ 
ability  literature  also  derive  the  availability  equations  fo**  two-unit 
standby  redundant  systems  with  the  multiple  repair  assumption.  Literature, 
is  scarce,  however,  for  two  unit  full-on  and  standby  Darallel  systems  with 
sinqle  repair  capability  H.e.  onlv  one  repair  crew  is  available  and  hence, 
failed  units  can  only  be  repaired  one  at  a  timel;  it  is  virtually  non¬ 
existent  with  respect  to  models  of  systems  made  uo  of  more  than  two 
redundant  units. 


This  paper  documents  the  derivation  of  availability  models  for  two 
unit  and  three  unit  full-on  and  standby  parallel  systems  with  sinqle  repair 
capabilities.  Also  included  are  derivations  for  three  unit  full-on  and 
standby  systems  with  sinqle  and  multiple  repair  capabilitv.jR|Jn  all  cases, 
each  unit  will  be  assumed  to  have  a  failure  rate  =  \  and  a  rpoair  rate  =  1 1. 

\ 

Case  1:  The  system  consists  of  two  units  both  fully  operational  whenever 
possible  and  has  sinqle  repair  caoability.  The  svstem  is  caoable  of 
performing  its  mission  when  onlv  one  unit  is  operational.  There  are  three 
possible  states  for  the  system  to  be  in  at  anv  given  time: 


State  2  - 

1  - 

0  - 

both  units  operating 

one  of  the  units  has  failed,  other  unit  operating 

both  units  have  failed. 

Using  a  Markovian  approach  and  letting  P2(tl,  P-|(t),  and  PQ(tl  represent 
the  probability  of  being  in  states  2,  1,  and  0  respectively  at  time  t 

yields: 

P2(t+At)  = 

P-|  ( t+At)  = 

P  (t+At)  = 

0 

P2(t)  H-2\At)  +  P 1  ( 1 1  jiAt 

P2(t)  (2\At)  +  P1f t)  (l-(\+n)Atl  +  PQ f 1 1  fiAt 

P1(t)-  (\Atl  +  PQ(tl  (1-uAt) 

Expanding  and  rearranging  terms  using  the  definition  of  a  derivative: 

p2<t)  * 

-2XP2(t)  +  [i  P,(t) 

p,(t> 

2\P2^tl  -  f\+nl  P^t)  +  |jiP0f t) 

p0u>  = 

XP^tl  -  |i  P0(tl 

Taking  Laplace  transforms  and  realizing  that: 

P2 ^ 0 )  =  1,  P -j ( 0)  =  0  =  PQ(0}  Assume  all  units  ooerative  at  t=0 


2 


f S+2X)  P?( S)  -  \i  P^S'  =  1 


2\  P2fS)  -  fS+X+tx)  P^Si  +  \i  PQ(S1  =  0 

X  P^S^  +  (S+[ii  PQfSl  =  0 


Solving  these  simultaneous  equations  for  PQ  results  in: 


p°'s'  3  — • - 

Sf S2  +  S(3X*2ti'  *  (2\2  +  2\H  +  |i2" 


The  denominator  must  first  be  factored  and  then  oartial  fraction  expansion 
accomplished  to  put  the  equation  in  a  form  to  find  the  inverse  Laplace 
transform. 


This  yields: 


P  (5^ 
o 


2\‘ 


where  S?  =  j[-(3X+?ji)  +  /\^+4Xn 


S^S-S-jUS-S^l  and  are  the  roots  of  S^+S(3\+?n)+(?\^+2Xn+n~) 


2\ 


S1S2S 


S^-S-jWS-S^ 


s?(s?-s1)(s-sp) 


3 


Findinq  the  inverse  transform  yields: 


P(t)  = 
o 


S,CS2-S,) 


S?(S2-S,) 


(S2-S,)  -  +  S,eV 


So  availability  then  is: 


A  =  l-PQ(tl 


S,S2iS2-S,)  -  2\2  [s2-S,-S2eSlt  *  S1eS?t 


sisats2-s1) 


SlV2V 


*  2\2  S?eSl^  - 


S1S2 


WVV 


2\n  +  n2  *  n2  l’s_eSIt  -  S,eS2 


2\~  +  2x4+11*  s1s?rs2-s1i 
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Case  2:  The  system  consists  of  two  units  with  one  unit  operated  and  one  7'n 
a  standby  state  and  has  single  repaid  capability.  When  a  failure  occurs, 
perfect  switching  is  assumed  to  connect  the  standby  unit  and  system 
operation  is  not  affected.  The  standby  unit  is  assumed  to  have  a  failure 
rate  equal  to  0,  the  operating  unit  is  assumed  to  have  a  failure  rate  equal 
to  \.  A  unit  under  repair  is  assumed  to  have  a  repair  rate  =  p.  The  three 
oossible  states  are  the  same  as  in  Case  1. 

Equations  of  state: 

P2( t+At)  =  P2(t)  (l-\At!  +  Pj(t)pAt 

P] ( t+At)  =  P2(t)  XAt  +  P^t)  (l-(\+plAtl  +  PQf t)pAt 

PQ(t+At)  =  P^tJXAt  +  P  (t)  M-nAtl 
P2ftl  =  -\P2ftl  +  pP^tl 

*•, ( t)  =  \P2(t)  -  fWP/ti  +  ||P  (t) 

P0(t)  =  XP^tl  -  n  Po(ti 
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Taking  Laplace  transforms: 


fS+\)  P2fSi  -  n  P}(S)  =  1 

p?(S)  -  rs+x+n)  P-,f s>  +  \l  pq> st  =  o 

P^S)  +  (S+fi)  PQ(S1  =  0 


Solution  for  PQf SX  after  solving  the  above  equation  set  simultaneously: 

A2 _  where  S-j,  S?  =  -(A+u)  + 

PQ(S1  =  S^S-S^fS-Sgl  These  are  the  roots  of  S^+Sf 2A+2n)+(A2+A^+p2) 


Having  the  denominator  factored  enables-  partial  fraction  expansion  and 
hence  finding  the  inverse  transform. 

The  Inverse  transform  of  the  above  results  in: 


Avai 1 abi lity : 


AC  t)  =  1-PJtl 


2  2 
V+u  +  r 


S2eV'  -  S,eV 


2  2 

r+xn+ji- 


SlVS2"V 


For  three  unit  systems,  there  a^e  4  Dossible  states: 


3  -  all  units  operable 

2  2  units  operable,  1  unit  failed 

1  -  1  unit  operable,  2  units  failed 

0  -  all  units  failed 


Case  3:  Three  different  three  unit  systems  will  be  investigated. 
Three  unit  full-on  parallel  operation  with  sinqle  repair. 

State  Equations: 


P3(t+At)  = 
P2(t+Atl  = 
P  •]  ( t+At)  = 
P  (t+At)  = 


P3(t)(l-3XAt)  +  P  ^ t^nAt 
P3(tl  3\At  +  P2(t)M-C2\+fxlAt)  +  P-jftlnAt 
P2ftl  2>At  +  P^tHl-f W&t'  +  PQC tlpAt 
P^tl  \At  +  P  (tUl-tiAtl 
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P3(tl 

=  -3\  P3(tl  +  p  P2(t' 

« 

P2(tl 

=  3\  P3ft)  -  (2\+p'  P2ft) 

+  P  P^t) 

p-,(t) 

=  21  P2(tl  -  ( X+pl  P^tl 

+  P  P0(t) 

P0(t) 

=  X  P,( t)  -  p  PQ(t) 

Taking  Laplace  Transforms: 


(S+3\)  P3(S )  -  p  P2(S)  *  1 

3\  P3(S)  -  ( S+2\+p)  P2r S)  +  p  P^Sl  =  0 

2\  p2(s)  -  fsn+p)  p^si  +  p  pQ(Si  =  o 

\  P^S)  -  (S+p)  P  (SI  =  0 

Simultaneous  Solution  of  the  above  for  P  (SI  results  in: 

o 

3 

P  (S)  =  -6\  where  S^,  S2,  S3  are  roots  of 

s(s-s1hs-s2hs-s3) 

S3+S2(5X+3p)  +  S(  ll\2+9^p+3p2i  +  (6\3+3^2+p3+6\V 
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To  find  values  for  S^,  $<>,  and  note  the  cubic  equation  is  of  the  form 


S3  +  pS2  +  qS  +  r  =  0 

6\  +  3fi 

in2  +  9\11  +  3  \12 
6\3  +  6\2|i  +  3\\i2  +  n3 
Substituting  (x  -  p/3)  for  S  reduces  the  equation  to: 
x3  +  ax  +  b  =  0 
Where:  a  =1/3  (3q-p') 

b  =  1/27  ( 2p3  -  9pq  +  27rl 

7  3 

Since  _b”+a_  <  0  there  are  3  distinct  real  roots. 

4  27 

Determine  the  value  of  the  9  in  the  expression  cos  9  = 
The  roots  x-j,  %2'  x3  will  have  the  values:  x-|  3  ^ 


Where:  p 


q 


r 


cos  £ 
3 
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‘2  *  2  "I!  +  f] 


x3  =  2  \f~JL  cos  f£  +  Ail 
3  l  3  3  . 


S2,  S3  are  then  calculated  from: 


X1  “f 


The  inverse  transform  for  P  is: 


P  rti  =  6\J  x 
0 


S,t  S?t  ^3"* 

fS1-S3^S1-S,US2-S3>-S2S3'S,-S3'e  +  ^VW61  -S  nS,  f  S  3 -S,!  e 


.  S1S2S3(SrS3WSrS2WS2-S3' 
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Avai labi litv: 


Aft'  -  1  -  P0(t' 

*  H3+3Xn3+6X3n  + 

ji3+3X^2+6X2h+6X3 

,r  sit  S-t  S  f 

6X3  S2S3(S2-S3'e  1  -  S1S3fSrS3'e  *  +  S^S^SyS^e  " 

SlS2VWfSl'S2U  W 

Case  4:  Three  unit  full-on  parallel  operation  with  multiple  reoair. 

State  equations: 

P3f  t+At'  =  P3^t'M-3\At'+P2ft'(iAt 
P  ft+At)  =  P3ft'  3XAt  +  Pgft'f  1-fPX+ji'At'  +  P^t'ZixAt 

P  i  ( t+At)  =  Pgf  t'  2XAt  +  P^t'n-fV^jx'At'  +  PQ't'3nAt 
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PQ(t+AO  = 

P^tl  \At  +  P0(tHl-3nAtl 

?3(t)  - 

-3\  P3ft)  +  |x  P2^tl 

P2(t' 

3\  P3(t)  -  ( 2\+p)  P2(t)  + 

Zp  P,(t' 

P](t)  = 

2\  P2(tl  -  f\+2n)  P^t)  + 

3n  P0ft) 

P0(t)  • 

X  P-,(t)  -  3u  PQftl 

Taking  Laplace  Transforms: 

( S+.3X)  P3(S1 

-  n  P2(S1  =  1 

3\  P3(S)  -  (S+2\+p)  P2(S)  +  2n  P}(s)  =  0 
2\  P2fS)  -  (S+\+2^  P ( S>  +  3(i  P^sl  =  0 
\  p  1  r s)  -  fs+3^)  PQfS)  =  o 


Simultaneous  Solution  of  the  above  for  P  (S): 

-6X3 

Po(S)  =  - 

S(S-S11(S-S2)(S-S31 

Where  S-j,  $p,  S3  are  roots  of  S3  +  S2  (6X+5ji!+SMH  \+\l)2)+6(  \+nl3 

S  i  a  + 

S->  =  -2(  X  +  n) 

S3  =  -3(\  +ji) 

The  Inverse  transform  is  the  same  as  for  Case  3  except  for  the  values  of 
Sp  Sp,  and  S3«  Availability  then  is: 

Aft!  =  1  -  P0(t) 

3  ?  ? 

=  \i  +3X»x  +3V \i  + 

(X+pO3 

3  T  S?t  S^t ' 

6\J  S2S3(S2-S3)e  -  S1S3fS1-S3)e  6  +  51S2(S-,-S2)e  3 

s1s2s3fs1-s3)fs1-s2>(s2-s3i 
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Derivations  for  the  standby  parallel  system  sinqle  multiple  ^eoa1' r  a^e 
found  in  the  same  manner  as  above. 


The  results  are: 

Case  5:  3  unit  standby  with  single  repair: 

A(t) 

3  2  2 

=  + 

3  2  .,2  ..3 

(1  +X|JL  +\  |i+\ 

,  s,t  s~t  s,t' 

\3[s2s3fS2-s3)e  -  S1S3^SrS3)e  +  S^fS^e 

S1S2S3  (SrS2USrS3VS2"S3^ 

Where  S-j,  S2,  S3  are  roots  of  S3+S‘'(3X+3|i)+S(3X^+4Xp+3p^l+(\3+\^p.+X|i“+p3l 

S-j  =  -(X  +  ul 

S2,S3  =  -(  X  +  |il  +\/2Xp 


Case  6:  Three  unit  standby  with  multiple  repair: 


Substituting  (x  -  p/3 ^  for  S  reduces  the  equation  to: 
x3  +  ax  +  b  =  0 
Where:  a  =  1/3  (3q-p*l 

b  =  1/27  l 2p3  -  9pq  +  27 r) 


2  3 

Since  Jl  +  JL  <  0  there  are  3  distinct  real  roots. 
4  27 


Determine  the  value  of  the  9  in  the  expression  cos  9  =  -b  f  2 


The  roots  x-j,  x^,  x^  will  have  the  values:  x-j  = 


2  cos  _9 
Y  3  3 


'*’2  &  C0$[V¥] 


x,  =  2  .  A&.  cos 


[fr] 


S-j,  S£,  S3  are  then  calculated  from: 


'I  3 


-f 


x  -  £■ 

x3  3 
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Documentation  for  the  derivation  of  all  the  availability  equations  tor 
redundant  systems  listed  in  Section  10  of  MIL-HDBK-338,  Electronic  Reli¬ 
ability  Design  Handbook,  could  not  be  found.  The  above  derivations  we^-e 
accomplished  to  correct  this  deficiency  and  to  check  the  accuracy  of  those 
equations. 
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